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ABSTRACT 


An  earlier  version  of  the  uncoupling  technique  referred  to  as  the 
Inertial-Damping  Collocation  Approximation  is  shown  to  be  improperly 
formulated,  and  a corrected,  improved  and  simplified  version  is  presented. 
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q - a"q  = -y”  + BQ^  + CQ^ 


•q  = DQ2  - 


Q = Ql  + Q2 


In  Eqs.  (l)-(3),  a dot  denotes  differentiation  with  respect  to  time  t,  a 
superscript  - 1 denotes  an  inverse,  and 

(i)  a and  y,  as  discussed  in  Ref.  [1],  are,  respectively,  known 
diagonal  damping  and  virtual  mass  matrices  corresponding  to 
appropriate  orthogonal  surface  expansion  functions. 

(ii)  a is  a matrix  intended  to  be  used  to  make  the  IDCA  go  over  to 
the  Curved  Wave  Approximation  (CWA)  for  short  times  (Ref.  [2]). 


1.  Earlier  Version  of  the  IDCA 


If,  for  simplicity  of  presentation,  it  is  assumed  that  the  excitation 

involves  no  incident  fluid  velocity  and  no  discontinuity  in  incident  pressure, 

the  Inertial-Damping  Collocation  Approximation  (IDCA)  presented  in  Refs.  [1] 

*) 

and  [2]  may  be  written  as 


(iii)  q is  the  matrix  of  (modified)  generalized  coordinates  for  the 
normal  shell  displacement  (corresponding  to  surface  expansion 
functions,  as  opposed  to  shell  modes). 

(iv)  and  are,  respectively,  inertial  and  damping  "components" 
of  the  generalized  force  Q corresponding  to  the  radiated  fluid 


pressure  on  the  shell  surface. 


White  StctloQ 
Buff  Stctlon 


It  should  be  noted  that  in  Eq.  (1),  and  in  what  follows,  the  sign  of 

the  virtual  mass  matrix  has  been  changed  to  produce  a positive  

definite  matrix.  L 
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(v)  B,  C and  D are  fitting  matrices,  to  be  determined  by  matching 
the  steady-state  response  of  a fluid  cavity,  having  the  same 
shape  as  the  shell,  at  selected  frequencies,  with  the  response 
obtained  using  acoustic  computer  codes. 


As  shown  in  Ref.  [2],  the  Doubly  Asymptotic  Approximation  (DAA)  may 


be  written  as 


-1*  -1 

q = -aQ-Y  Q 


which,  for  short  times  (or  high  frequencies),  becomes  the  Plane  Wave 


Approximation  (PWA) 


Q = -a  q 


and,  for  long  times  (or  low  frequencies) , goes  to  the  Virtual  Mass 


Approximation  (VMA) 


as  it  should. 


Q = -Y  q 


If  B,  C,  D and  a are  set  equal  to  zero,  Eqs.  (l)-(3)  lead  to 


Q = -a  q - Y q 


which  does  not  go  to 


although  using 


Furthermore, 


a - -C  8 

in  Eq.  (1),  as  in  Ref.  [2],  where  8 is  a symmetric  curvature  matrix  defined 
in  Ref.  [2],  leads  to  the  CWA  when  B and  D are  set  equal  to  zero,  an 
inertial  component  has  been  used  to  satisfy  short-time  damping  behavior. 

This  would  cause  a very  poor  match  at  high  frequency  in  the  relation 
between  the  generalized  force  and  the  inertial  component  of  the  acceleration 
as  a function  of  frequency. 


J 
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2.  Corrected  Version  of  the  IDCA 


Instead  of  Eqs.  (l)-(3),  an  appropriate  form  is 

.-1. 


- BX  "q^  = - Y Q + B Q 


q‘2  - Deq2  = DQ  - a Q 


^ * "2 


in  which  X and  e are  matrices  introduced  to  permit  the  IDCA  to  approach 
the  correct  asymptotic  limits.  If  B and  D are  set  equal  to  zero,  Eqs,  (9)- 
(11)  go  to  the  DAA  of  Eq.  (4)  as  they  should. 

3.  Determination  of  B and  D 

On  the  surface  of  a fluid  cavity  having  the  same  shape  as  the  shell, 
apply  the  pressure 

/s 

p^(s,  t)  = t|»^(s)e^  ^i*^ 

In  Eq.  (12),  s denotes  position  on  the  shell  surface,  is  a constant, 
i-/^,  is  an  appropriately  chosen  fitting  frequency,  and  ip^(s)  is 
one  of  the  orthogonal  surface  expansion  functions  which  satisfy 


■ “1 


where  A is  the  surface  area  of  the  shell  and  6^^  is  the  Kronecker  delta. 
The  generalized  forces  corresponding  to  the  pressure  of  Eq.  (12)  are 


«j>i 


p^(s,  t)  'l'j(s)  dA 


Using  Eqs.  (12)  and  (13),  Eq.  (14)  becomes 


(9) 

(10) 

(11) 

■j 

i 

j 

i 

I 

4 

•1 

(12)  I 


(13) 

(14) 


(15) 
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which  yields 


in  t 

(Ql)i  = = Q^(t) 


(Qj)j^  = 0,  j ,t  i 


On  the  surface  of  the  shell,  the  normal  acceleration  resulting  from 
the  excitation  of  Eq.  (12)  may  be  expressed  as 

w^(s,  t)  - Q^(t)A^(s)  - Q^(t)A^(s) 

Since  the  normal  acceleration  may  also  be  written  in  the  form 


= y iJ)  (q  ) 
k 


it  follows  that 


■ 

I «k)l  Vj  Vj  ■iA  - Qi 

J ) 


In  view  of  Eq.  (13),  Eq.  (19)  becomes 


(qj)i  * ■ ^i^'^2j^l 


where 


(’'lj>l  * t J 

^ A 


J A 


In  order  to  avoid  oscillations  in  the  relation  between  generalized 
forces  and  acceleration  obtained  from  Eqs.  (9)-(ll)  as  a function  of 
frequency,  as  discussed  in  Ref.  [1],  the  (Wj^j)j^  and  terms  (j  / i) 

will  be  neglected  compared,  respectively,  to  the  and  (^2^^)^^  terms. 


T 
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This  assumption  neglects  coupling  between  surface  expansion  functions  at 
all  frequencies  and  allows  Eq.  (20)  to  be  simplified  to 


^1  “ "^l^^ll  ^1^^21 


In  which 


RllC^l)  - . *^21^^1^  * ^^^^^21^1 


The  first  term  on  the  right-hand  side  of  Eq.  (22), 


qii  = 

Is  the  real  (Inertial)  part  of  the  second  derivative  of  the  generalized 
coordinate.  In  phase  with  the  applied  pressure,  and  the  second  term. 


^21  ‘ V^21 


Is  the  Imaginary  (damping)  component,  90°  out  of  phase.  It  now  remains 
to  match,  at  selected  frequencies,  the  results  yielded  by  Eqs.  (9)- (11) 
with  those  yielded  by  Eqs.  (24)  and  (25). 

Consider  simple-harmonic  motion,  at  frequency  Involving  the 

1th  surface  expansion  function  only.  If  the  matrix  X Is  assumed  to  be 
diagonal,  Eq.  (9)  may  be  satisfied  by  a diagonal  matrix  B with  main 
diagonal  terms  such  that 

^11  " Vl  *^11  “ “®1^1  * ®1^1 

In  which  e^,  the  1th  main  diagonal  term  of  the  diagonal  matrix  y^.  Is 


given  by 


'i  ■ 


where  y^  Is  the  1th  main  diagonal  term  of  the  virtual  mass  matrix.  Since 
a steady-state  excitation  Is  being  considered,  Eq.  (26)  may  be  written 


- 6 - 


in  the  form 


’ll  - + “l\'‘ll> 


In  order  for  Eqs.  (24)  and  (28)  to  yield  identical  results,  at  frequency 


«ii. 


-QlC.l  + 


or 


1 - e R ((2  ) 

B ^ H B ((2  ) 

1 fn  /r\  \ •v  in^  i -Li 


t^i^^ii>  - 


Consider  now  simple-harmonic  motion,  at  frequency  ^22^,  involving 
the  ith  surface  expansion  function  only.  If  the  matrix  e is  assumed 
to  be  diagonal,  Eq.  (10)  may  be  satisfied  by  a diagonal  matrix  D with 


main  diagonal  terms  such  that 


^2i  ~ °i^i^2i  °i^i  ^i^i 


-1 


in  which  h^,  the  ith  main  diagonal  term  of  the  diagonal  matrix  ot  , 


is  given  by 


hi  = 1/a^ 


where  0(^  is  the  ith  main  diagonal  term  of  the  damping  matrix.  For  the 


steady-state  excitation  being  considered,  Eq.  (30)  may  be  written  as 


’21  ■ A*"!  ^ * Vl'4l> 


In  order  for  Eqs.  (25)  and  (32)  to  yield  identical  results,  at  frequency 


^2i’ 


-Qi(h^  + D^/J22^)/(1  + - -Qi/R2i^“2i^ 


or 


[1  - h.R-.((2,.)]n,, 

D ■ * ■* “ »D  ((2  ) 

i ^21^^21^  ■ ^1  “ i 2i^ 


(28) 


(29) 


(30) 


(31) 


(32) 


(33) 
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4.  Positive  Inertia  and  Damping  - Allowable  Fitting  Frequencies 

For  the  steady-state  excitation,  at  frequency  of  the  ith  surface 
expansion  function  only,  Eqs.  (9)-(ll)  may  be  written  as 
[e^  + [h^  + 

^ [1  + ^ [1  + 1 


(34) 


in  which  and  given  by  Eqs.  (29)  and  (33),  respectively. 

If  and  do  not  vanish,  Eq.  (34)  reduces  to  the  PWA  as  S)  -*■  «>: 


Him  ^ 


and  reduces  to  the  VMA  as  fi  -►  0: 


Jlim 
^i  " fW) 


1 


-i  Qi 


If  Eq.  (34)  is  to  produce  positive  inertia  and  damping  at  all  frequencies 
(2,  the  coefficients  of  and  must  be  negative.  Thus,  if  X^  and  are 
positive,  positive  inertia  and  damping  are  ensured  by  requiring  B^((2j^^), 
^i^^2i^  > 0.  From  Eqs.  (29)  and  (33)  it  then  follows  that 


L - eiRii(^ii) 
^li^^li^  “ ^i 


> 0 


1 ^i^21^Si^ 

^21^Si^  - 


> 0 


(35) 

(36) 


are,  respectively,  sufficient  conditions  for  positive  inertia  and  damping. 

The  inequality  of  Eq.  (35)  will  be  satisfied  if  the  fitting  frequency 
is  such  that  Rj^j^(^2,i^  between  and  X^,  i.e., 

0 < < Yj 


(37) 


n ' I 
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Similarly,  the  inequality  of  Eq.  (36)  will  be  satisfied  if  the  fitting 
frequency  is  such  that  *^2i^^2i^  lies  between  and  e^,  i.e., 

0 < < R2i^^2i^  ^ “i 

Since  X is  the  high  frequency  limit  of  the  inertial  part  of  Eq.  (34) 

1 1 

•J 

and  is  the  low  frequency  limit  of  the  damping  part  of  Eq.  (34)  (both 
of  which  limits  vanish  in  an  exact  solution) , it  is  recommended  that  and 
be  arbitrarily  set  to  small  fractions  of  end  , respectively.  The 
fractions  may  be  determined  by  comparing  the  frequency  behavior  of  the  inertial  ! 

and  damping  parts  of  Eq.  (34)  with  those  of  Eq.  (22).  Least  squares  fits 
of  the  type  proposed  in  Ref.  [3]  may  also  be  used  to  detennine  the  X^  and 
£^.  Alternatively,  experience  may  Indicate  that  the  matrices  X and  e should 

be  chosen  such  that  known  steady-state  results  from  acoustic  computer  codes  ! 

I 

are  matched  at  low  and  high  frequencies,  respectively,  while  still  maintaining  j 

the  inequalities  of  Eqs.  (37)  and  (38).  In  any  event,  the  satisfaction  of 
the  inequalities  of  Eqs.  (37)  and  (38)  ensures  positive  inertia  and  damping. 

If  the  inequalities  of  Eqs.  (37)  and  (38)  are  satisfied,  it  may  be 
readily  shown  that  in  a steady-state  excitation  at  any  intermediate  frequency 
0(0  < 0 < ”) , the  inertia  and  damping  yielded  by  Eq.  (34)  [or  equivalently 
Eqs.  (9)-(ll)]  are  smaller  than  the  corresponding  quantities  produced  by 
the  VMA  and  the  PWA.  A closer  reproduction  of  the  exact  curve  relating 
Inertial  generalized  force  and  driving  frequency  might  be  obtained  by 
introducing  an  additional  fitting  matrix  in  Eq.  (9),  enabling  an  additional 


point  of  the  approximate  curve  to  match  that  of  the  exact. 
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